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摘 要 研究 了 一 道 工序 智能 加 工 系统 轨道 式 自动 引导 车 (RGV) 调度 问题 . 该 问题 为 
2018 年 全 国 大 学 生 数 学 建 模 竞赛 B 题 的 一 部 分 . 系统 由 一 辆 轨道 式 自动 引导 车 和 若干 台 计 
算 机 数控 机 床 (CNC) 等 部 件 组 成 , RGV 操控 多 台 CNC 完成 多 个 物料 加 工 , RGV 调度 方案 决定 
了 系统 的 效率 . 以 RGYV 的 移动 路 径 为 决策 变量 ,以 RGV 在 CNC 上 的 操作 结束 时 刻 为 时 间 
节点 ， 以 物料 加 工 剩余 时 间 为 状态 变量 , 给 出 了 问题 的 数学 模型 , 但 模型 中 的 部 分 参数 以 决 
策 变量 为 下 标 . 通过 定义 新 的 变量 和 约束 ,将 模型 修改 为 不 含 变量 下 标 和 分 段 函 数 的 非 线性 
混合 整数 规划 模型 . 最 后 给 出 了 算 例 ， 说 明了 模型 的 正确 性 和 可 操作 性 . 
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with One Process 
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Abstract: Scheduling problem of intelligent processing system is studied. This problem is a part 
of Problem B of 2018 China Undergraduate Mathematical Contest in Modeling. The System 
consists of a Rail Guide Vehicle (RGV), several Computer Number Controllers (CNC) and other 
components. RGV manages multiple CNCs to finish multiple units of material. RGV scheduling 
scheme determines the efficiency of the system. Taking RGV's moving path as decision variable, 
RGV's operation ending time on CNCS as time nodes, and material processing remaining time as 
state variables, the mathematical model of the problem is developed. However, the subscripts of 
some parameters are decision variables in this model. By defining new variables and constraints, 
the model is modified to exclude the decision-variable subscripts and piecewise functions, and the 
model is transformed into a nonlinear mixed integer programming model. Finally, a numerical 
example is given, which illustrates the correctness and operability of the model. 
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给 定 的 智能 加 工 系 统 由 一 辆 轨道 式 自动 引导 车 (Rail Guide Vehicle, RGV) 和 若干 台 计算 机 
数控 机 床 ( Computer Number Controller CNC) 等 部 件 组 成 ， RGYV 操控 多 台 CNC 完成 多 个 物 
料 加 工 ， RGV 调度 方案 决定 了 系统 的 效率 . 该 问题 为 2018 年 全 国 大 学 生 数 学 建 模 竞赛 赛 题 
[1], 文献 [2 指出 该 问题 是 一 个 具有 很 强 的 实用 性 和 操作 性 的 问题 

系统 可 描述 为 : RGV 是 一 种 无 人 芍 驶 、 能 在 固定 轨道 上 自由 运行 的 智能 车 , 根据 指令 
能 自动 控制 移动 方向 和 距离 ,并 自 带 一 个 机 械 手 臂 , RGV 的 机 械 手 臂 前 端 有 2 个 手 爪 ， 通 过 
旋转 可 以 先后 各 抓 取 1 个 物料 , 完成 上 下 料 和 清洗 熟 料 作 业 ,， 同 一 时 间 只 能 执行 移动 、 停 目 
等 待 、 上 下 料 和 清洗 作业 中 的 一 项 ; 每 台 CNC 安装 同样 的 刀具 , 物料 可 以 在 任意 一 台 CNC 
上 加 工 完成 ， 物料 加 工具 需 一 道 工 序 ; 传送 带 为 CNC 输送 生 料 (未 加 工 的 物料 )， 将 成 料 


上 
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《加工 并 清洗 完成 的 物料 ) 送出 系统 . 只 考虑 RGV 移动 、 上 料 、 下 料 和 清洗 熟 料 的 时 间 对 
系统 效率 的 影响 . 关于 系统 的 更 详细 的 介绍 参见 文献 [1]. 

文献 [2] 先 分 析 一 个 周期 RGV 运行 路 径 的 耗 时 ， 并 给 物料 编号 ， 以 某 个 物料 是 否 指派 到 
某 个 CNC 为 0-1 决策 变量 ,基于 耗 时 最 小 路 径 建立 了 非 线性 混合 整数 规划 模型 ， 设 计算 法 求 
解 模型 .本文 给 出 了 与 文献 [2] 不 一 样 的 建 模 方法 . 文献 [3] 以 仿真 模拟 的 方法 研究 了 该 问题 . 
文献 [4] 用 动态 规划 方法 研究 了 一 个 周期 RGV 运行 路 径 最 短 的 问题 . 
1 建 模 
1.1 符号 和 分 析 

一 个 物料 的 加 工 过 程 为 : 上 料 , CNC 加 工 物 料 ， 下 料 ， 最 后 是 清洗 , 其 中 上 料 、 下 料 和 清 
洗 都 需要 RGV 来 操作 . 下 料 后 , RGV 可 以 先 清洗 物料 ,也 可 以 通过 男 一 手 爪 ， 先 完 成 男 一 物 
料 的 上 料 操 作 ， 再 清洗 物料 . 不 难 发 现 先 完成 另 一 物料 的 上 料 更 节约 时 间 . 一 辆 RGV 要 同时 
操控 m 台 CNC, RGV 从 CNCi 移动 到 CNC j 需要 的 时 间 为 di (distance). 每 个 物料 的 加 工时 间 


一 


为 p (process)， 清洗 时 间 为 w (wash)， 当 p<w 时 只 需要 一 台 CNC 工 作 , 所 以 假设 p>w. 
CNCi 的 上 料 时 间 为 a (arrive)， 下 料 时 间 为 1 (eave). 其 它 操作 时 间 忽 略 不 计 . 


在 没 被 物料 占用 的 CNC 上 , RGYV 的 操作 为 上 料 一 次 ; 在 被 占用 的 CNC 上 ,RGYV 的 操作 
为 先 下 料 , 然后 上 料 , 最 后 清洗 , 将 在 同一 CNC 上 连续 的 下 料 、 上 料 和 清洗 也 记 为 一 个 操作 ; 
xz=0 表 示 给 空闲 的 CNC 上 料 ，z=1 表 示 给 加 工 完 的 CNC y, 下 料 、 上 料 , 然后 清洗 物料 ;n 
个 物料 对 应 n 个 上 料 操作 . 不 再 取 生 料 的 时 刻 之 后 , CNC 上 仍 有 没完 成 加 工 的 物料 , RGV 到 
达 各 个 CNC 时 , 需要 执行 的 操作 为 : 下 料 ， 然 后 清洗 熟 料 . 调度 RGV 的 顺序 可 以 按照 先 加 
工 则 先 下 料 的 顺序 执行 , RGV 完成 一 个 CNC 上 的 操作 ， 就 开始 往 下 一 个 CNC 移动 . 为 了 模 
型 的 简洁 , 我 们 在 模型 中 不 考虑 不 再 取 生 料 的 时 刻 之 后 的 加 工 过 程 ， 以 RGV 从 开始 工作 到 
第 n 个 操作 结束 时 刻 的 时 间 段 最 短 为 目标 . 

RGYV 的 初始 位 置 为 y,， 按 时 间 先 后 顺序 的 n 个 操作 对 应 RGV 的 n 个 位 置 , 记 为 RGV 


的 移动 路 径 y=(y,y,,…,y,). 称 RGV 在 CNC y, ,的 操作 zz 结束 时 刻 到 在 CNC y 的 操作 z， 


结束 时 刻 为 阶段 i. 阶段 i, RGV 先 从 CNC y ,移动 到 CNC y,, 然后 在 CNC y, 上 对 物料 执行 


操作 . 在 RGV 到 达 时 ,如 果 CNC y, 加 工 物 料 还 没完 成 , RGV 还 需 等 待 . 
阶段 i 开始 时 刻 ( 即 阶段 i-1 结束 时 刻 )CNC 的 状态 向 量 记 为 


hi > (hi11 hi A hi ), 


用 及 =-1 表 示 阶 段 i 结束 时 刻 CNC j 上 没有 物料 , 对 应 RGV 的 操作 为 =0; h, >0 时 , 阶 


段 i 结束 时 刻 CNC j 上 的 物料 加 工 还 需 h, 单 位 时 间 , 对 应 RGV 的 操作 为 z, =1, 表示 有 加 工 


完成 的 物料 在 等 待 . 初始 时 刻 为 0, 初始 状态 为 h =(-l, 一 L…, 一 ]). 


阶段 i, 给 定 状 态 h_ ,和 RGYV 位 置 y ,决策 为 确定 RGV 的 下 一 位 置 wy 和 操作 z,. 因 决 
2 


态 hy 


1.2 阶段 时 长 

RGYV 到 达 CNC 时 ， 
果 CNC 空闲 ，R 
的 情 


攻 见 表 1 的 第 三 


CNC 状态 


1.3 ”状态 转移 公式 


对 于 CNCy,: 若 有 h 


i—1,y; 


物料 剩余 加 工时 间 为 p， 即 为. =p; 车，, >0， 则 操作 为 下 料 、 上 料 、 再 清 


RGYV 只 给 CNC 上 料 ， 见 表 1 的 第 二 行 。 
、 四 行 。 


确定 ，y 为 唯一 独立 的 决策 


根据 CNC 的 状态 执行 相应 的 操作 。 例 如 ，RGYV 到 达 CNC 时 ， 如 


表 1 CNC 状态 与 RGV 的 操作 
RGYV 操作 


空 朵 


| 
剩余 加 工时 间 大 于 零 Ee 
剩余 加 工时 间 等 于 零 有 


由 表 1, 不 难得 出 阶段 i 时 长 x 表示 为 


d, 1 六 十 0 ， 和 元 一 

A 三 4 ， L, Q， ,0<h, <d, a 
hi, L 4 四 hy, >d, 1 
= 一 1]， 即 CNC y, 上 没有 物料 , 操作 为 上 料 ， 


RGYV 到 达 CNC 时 ，CNC 上 有 物料 


i 


(1) 


上 料 结束 后 , CNC y 的 


后 , CNC y 的 物料 剩余 加 工时 间 为 pw， 即 及 =p 一 w. 


洗 , 操作 结束 


对 于 CNCJj ，7j zy: 若 有 ,j=-1, 则 hh =-1, 即 CNC 上 仍 没有 物料 ; 若 加 ,>0， 则 
hh; =max{0,h ,一 zw}， 即 CNC 上 的 物料 加 工 又 过 去 了 ,单位 时 间 . 阶段 i 结束 时 刻 的 状态 
hh 表示 为 

pb, hi;=—1, j= y,, 
PpP—wW, hi; ) 过 0,7 = yi 
h, —], 万 二 =—l, jz y,, ® 
max{0,h 1 ;—A)}, hi;>0,jzy 
1.4 模型 

以 时 间 长 度 最 小 化 为 目标 得 模型 1. 

模型 1 


201907.00015v1 


chinaXiv 


2 模型 转化 
因 含 变量 下 标 和 分 段 函 


定义 = (NX, XN, ), 二 


| 
d, ,+ a,, 1Ly 三 一 上 
A;=1d,,, +l, +a, ,Ozh i, <d,,,, 
hy, tl +a, +w, hi, >d,,,; 
pb, hj; = -1, = yY,, 
Ss.t. = P , 0 > yi 
-1, Di) = 一 
max {0,h 1 ;—NA)}, hj;>0,7zy; 
hh = (一 一 一 ,0 = 二 
y {1,2,.…,m}; 
i=1,2,.…,7;]=1,2,.…,m. 


数 的 模型 在 很 多 软件 中 不 能 直接 求解 ， 下 面 修改 模型 1 对 (]) 式 , 
(Go ， 满足 


> =L2 e{0,1}. 


1i=1 


Xj =1 表 示 RGV 在 第 i 阶段 的 操作 在 CNC j 上, 决策 由 y 改 为 X. 定义 e=(ei,es,en), 满足 


对 应 (1) 式 的 三 个 部 分 . 约束 


使 得 : 当 iu =1、 Xi, =1、 


示 足 够 大 的 整数 ,约束 


使 得 : 当 x ,=1、 罗 =1、 


使 得 : 当 z vs =1、 六 =1、 


7 = > [ej x; 1 X, (d,, 本 a,) 十 ep (CQ 


Be 


i—l,u™ iv “il i—l,u” iv i—l,v 


六， = -1， 一 定 有 e, =1. 用 表示 可 以 任意 小 的 正 数 ，M 表 


Mx, Xi Ci2 全 Xi_Lu Xiy (hi, 十 E)(d hs, 二 2) 


1 一 1 iv i 


万 1， (d hi,) 和 0 ( 即 0 全 hy < d, )， 一 定 有 Ci 二 1. 约束 


MXx, , Xx, es > xX, ,xX, (2 -d, 


1 一 1 iv Ti3 i—l,u” iv tv 


hs, —d,, > 0 入 一 定 有 Ci3 =1. (1) 式 修改 为 
+l,+a,+w) 


La 


t+E3Xi Lu Xiy (hi, 中 L, 让 人 中 w)]. (3) 


对 (2) 式 , 定义 w = (onsop)， 满足 


y 


omd+op =1,0;, € {0,1}. 


il 


对 应 hh ,= 一 1 和 为 ， ;>0 两 种 情况 ， 类 似 的 分 析 可 将 (2) 式 修改 为 
h,=onxyp+oXy(p—w)+on(l—X)(-D+o;(1—X%)max{0,h,, ;— x} (4) 
由 (3) 和 (4), 模型 1 修改 为 模型 2. 
模型 2 
和 
2 =1,x, e{0,1}; 
> =1,e Ci e {0,1}; 
Xi LuXivEil = Xi- Lu Xv hi lv? 
Mx,_ Lu XvEi2 > lu 如 (有 -Ly 0 和 h_ lv + €); 
Max, LuXivEi3 > Lu iy Xi,(h, -lv d,, —€); 
Sl Ti [eax lu Xd +a y) + CaN lu Xd +h, +a, ,+ W)+ EisXi_ lu Xi,(h, -ly +h, +a, + w)]; 
u=l] v=] 
e+e;s =1,0; {0,1}; 
Mo >—h_ 
Mo >h_ 于 
h, j = Op1Xy P+oXi(p W) 1+ oll Xj)( 1) + ll x;)max {0, hh 1 一 Z| 
hh = (一 上 一 一]D,2xo = (1,0,…, 0); 
i=1,.…,n;u,y, J =1],.…,m;s = 1,...,3;1 二 1,……,2. 
3 求解 


对 于 模型 1, 大 决 策 y 是 确定 的 , 则 可 以 从 第 一 阶段 开始 ， 由 式 (1) 和 (2) 循 环 求 x, 和 的 
目标 函数 值 ， 也 能 很 快 地 由 模型 1 求 出 在 “每 个 周期 RGV 移动 的 路 径 回 


值 ， 


定 ” 假 设 下 的 最 优 解 . 假 t 
第 三 组 数据 进行 计算 , 380 个 阶段 结束 时 刻 为 7.9236 小 时 , 文献 [2] 给 出 的 是 8 小 时 内 完成 393 
个 物料 的 加 工 ， 


能 很 快 地 求 出 


\ 若 


通过 


没 每 个 周期 RGV 移动 的 路 径 固定 ， 


过 Matlab 对 文献 [1] 的 算 例 的 


因为 我 们 的 模型 不 考虑 最 后 几 个 物料 的 加 工 过 程 ，380 个 阶段 结束 时 刻 CNC 


上 还 有 8 个 物料 在 加 工 , 所 以 我 们 的 结果 与 文献 [2] 给 出 的 结果 基本 一 致 . 由 zx, 入 的 值 不 


> 


较 多 非 线 
4 结语 


在 得 出 上 下 料 等 操作 的 开始 时 刻 , 
对 模型 2, 通过 Lingo 计 入 
生 约 束 ， 计 算 速 度 较 慢 ， 对 模型 2 设计 更 快 的 算法 是 需要 进 


了 只 有 2 到 5 个 物料 的 算 例 ,与 模型 1 的 结果 相符 . 因为 模型 2 含 


步 研究 的 课题 


时 间 节 点 和 状态 变量 的 选取 对 模型 1 的 建立 非常 关键 , 模型 1 不 需 先 考虑 一 个 周期 


RGY 的 耗 时 ， 


十 六 > 


它 CNC 的 状态 为 当前 状态 ， 和 


合 中 去 掉 ， 


其 重新 加 入 CNC 集合 , 并 将 其 状态 定义 为 空间 ， 


忽略 最 后 几 个 工件 的 做 法 也 使 得 模型 更 为 简洁 . CNC 故障 时 , 将 其 从 CNC 集 
出 现 故 障 的 CNC 修复 后 , 将 
它 CNC 的 状态 为 当前 状态 , 然后 重新 求 
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解 模型 . 加 工时 间 很 短 时 , 最 优 解 会 有 CNC 一 直 空闲 . 也 可 基于 模型 1 考虑 多 道 工序 的 问 
题 . 
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